Unconventional spin freezing and fluctuations 
in the frustrated antiferromagnet NiGa2S4 
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Muon spin rotation (/iSR) experiments reveal unconventional spin freezing and dynamics in the 
two-dimensional (2D) triangular lattice antiferromagnet NiGa2S4. Long-lived disordered Ni-spin 
freezing (correlation time > 10~® s at 2 K) sets in below T/ = 8.5 ± 0.5 K with a mean-field- 
like temperature dependence. The observed exponential temperature dependence of the muon spin 
relaxation above Tf is strong evidence for 2D critical spin fluctuations. Slow Ni spin fluctuations 
coexist with quasistatic magnetism at low temperatures but are rapidly suppressed for fields > 
10 mT, in marked contrast with the field- independent specific heat. The ^SR and bulk susceptibility 
data indicate a well-defined 2D phase transition at Tj, below which NiGa2S4 is neither a conventional 
magnet nor a singlet spin liquid. 

PACS numbers: 75.20.Hr, 75.40.Cx, 75.50.Ec, 76.75.+i 



The layered chalcogenide NiGa2S4 is an example of a 
quasi-two-dimensional (2D) antiferromagnet (AFM) with 
a nearly perfect triangular lattice, and is a candidate 
for a geometrically-frustrated spin liquid.— Its crystal 
structure consists of triangular Ni planes well separated 
by GaS polyhedra, which assure the quasi-2D nature 
of the Ni magnetism. The Ni^+ ^2g^'g configuration is 
that of a S* = 1 Heisenberg magnet, consistent with 
the nearly isotopic susceptibility. Neutron scattering 
experiments^''^ revealed the development of quasistatic 
(correlation time Tc > 3 x 10^^° s) short-range incom- 
mensurate Ni spin correlations below ~20 K, well below 
the paramagnetic Curie- Weiss temperature \9w\ ~ 80 K, 
with very weak correlations between Ni planes. The 
magnetic specific heat Cm(T) exhibits a broad peaki at 
~10 K, and the magnetic susceptibility shows a sharp 
kink at 8.5 K suggesting a transition at this tempera- 
ture. Evidence for spin freezing below ^10 K was ob- 
tained from ^^Ga NQR experiments'^ and substantiated 
by recent muon spin rotation (/zSR) datai^ 

In the temperature range 0.35-4 K Cm{T) oc T^, in- 
dicative of gapless linearly-dispersive low- lying modes. 
Magnons in an ungapped 2D AFM would yield this re- 
sult but would require long-range spin order, in appar- 
ent confiict with the neutron scattering results. At all 
temperatures Cm (T) is independent of applied magnetic 
field up to 7 T, as also observed in the 2D kagome lat- 
tice AFM SrCr9pGai2_9pOi9 (SCGO) (Ref.H). This is 
certainly not expected for magnons or any other simple 
cooperative spin excitation, and suggests the possibility 
of a singlet or singlet-like spin liquid. A number of theo- 
retical scenarios have been proposed for NiGa2S4 (Refs. ^ 
and0), with magnetic properties that depend sensitively 
on details of the assumed model as is common in frus- 
trated systems. It is important, therefore, to characterize 



the properties of this compound as fully as possible. 

This Letter reports a detailed longitudinal-field fiSR 
(LF-^SR) (Ref. 8) study of NiGa2S4, using the local na- 
ture of the muon spin probe to provide unique informa- 
tion on the static and dynamic behavior of Ni spins at 
low temperatures. Although preliminary /iSR results^ 
seemed to be inconsistent with spin freezing, the present 
data and those of Ref. H give clear evidence for a dis- 
ordered quasistatioi^ Ni-spin configuration. We find a 
lower bound on Tc more than three orders of magnitude 
longer than the neutron scattering value. This spin freez- 
ing sets in abruptly and exhibits disordered mean-field- 
like behavior below a freezing temperature Tf w 8.5 K. 
A kink is observed in the bulk susceptibility at the same 
temperature; together these results suggest a well-defined 
2D phase transition. Dynamic muon spin relaxation is 
spatially inhomogeneous at all temperatures, exhibits an 
exponential temperature dependence for T > Tf as ex- 
pected for 2D Heisenberg critical fluctuations, and re- 
mains strong down to ~2 K for small applied fields. Our 
most unexpected finding is a significant suppression of 
this low-temperature relaxation by fields 10 mT, which 
is difficult to reconcile with the field independence of the 
specific heat. 

Time-differential LF-/iSR measurements were carried 
out at the M20 beam hue at TRIUMF, Vancouver, 
Canada, on a homogeneous and stoichiometric polycrys- 
talhne powder and a single crystal of NiGa2S4 (Ref. [l|), 
the latter grown using a chemical vapor transport 
method. Data were taken in a longitudinal field fioH^ > 
2 mT (Hl II [001] for the single crystal), to decouple the 
muon spins from nuclear dipolar fields^^ and to study 
the field dependence. Care was taken to minimize the 
spectrometer "dead time" between a muon stop and the 
earliest detection of the decay positron. This was nec- 
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essary since the early-time relaxation in NiGa2S4 is very 
rapid; it was not observed in earlier /xSR experiments at 
the KEK pulsed muon facility^ due to the relatively long 
dead time inherent in pulsed /iSR. 

Figure [T] shows representative LF-/iSR asymmetry de- 
cay data^ in polycrystalline NiGa2S4. The asymme- 
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FIG. 1: (color online) Representative (a) early- and (b) late- 
time LF-/.iSR asymmetry data in polycrystalline NiGa2S4, ap- 
plied longitudinal field hoHl ~ 2 mT. Curves: fits to Eqs. ([l]- 
13). 

try A{t) is proportional to the muon spin polariza- 
tion P{t). In agreement with previous results,^ the 
early-time data at low temperatures exhibit strongly 
damped oscillations, indicative of a quasistatic compo- 
nent (Bloc) of the local field at the muon site. The os- 
cillations are rapidly damped [Fig. [Ija)] due to inhomo- 
geneity in (Bioc), leaving a second signal that relaxes 
slowly [Fig. mb)]. This slow signal is due to the muon 
spin component parallel to (Bioc) (Kubo-Toyabe behav- 
ior), and its relaxation is due to the fluctuating compo- 
nent SBiocit) = Bioc(0 - (Bloc) (Ref.ini). 

It should be noted that this two-component Kubo- 
Toyabe structure arises when the statistical properties 
of Bloc [inhomogeneity in (Bioc), fluctuation statistics 
associated with S'Biodt)] are the same at all muon sites, 
i.e., when the system is macroscopically homogeneous. 
Thus the two components should not be associated with 
separate volume fractions of a multiphase sample, and 
are not by themselves evidence for multiple phases in 
NiGa2S4. We discuss below independent evidence for in- 
homogeneous spin dynamics obtained from the form of 
the slow-signal relaxation. 



The oscillations and Kubo-Toyabe behavior are strong 
evidence for a quasistatic Ni-spin configuration in 
NiGa2S4. As Tf = 8.5 ± 0.5 K is approached from be- 
low, the oscillation frequency decreases rapidly and the 
two-component behavior is lost. For T > Tf there is only 
a single signal component, with a rapidly-increasing rate 
as T ^ Tf as expected from critical slowing of Ni-spin 
fluctuations. 

For T < Tf the asymmetry data were fit to the form 

A{t) = ArPrit) + A,P,{t) , Pr{0) = P,(0) 1 , (1) 

where the the first and second terms represent the rapid 
and slow components, respectively ("T < Tf fits"). The 
total asymmetry A = Ar + A^ at t = was assumed 
independent of temperature and applied field. A damped 
Bessel function 



Pr{t) = exp(-Art) Jo(t^^t) , 



(2) 



indicative of an incommensurate spin structure)^ was 
found to fit the early-time data better than other can- 
didates such as a damped sinusoid exp(— A^t) cos(a;^t) 
or simple Kubo-Toyabe functions . ^^'^"^ The spectrometer 
dead time of '^8 ns makes the choice of fit function un- 
certain, but the fit value of the oscillation frequency lu^ 
does not depend strongly on this choice. For the Bessel- 
function fit w^/27r = 34 ± 2 MHz, (Bioc) = ^^m/Tm = 
251 ± 15 mT, where 7^ is the muon gyromagnetic ra- 
tio. This value is qualitatively consistent with Ni dipolar 
fields at candidate muon stopping sites in the GaS layers. 
The ratio Xr/ujfj. ~ 0.2 is a measure of the relative spread 
in (Bloc) due to disorder. 

The late-time asymmetry data below Tf and the en- 
tire asymmetry function above Tf {Ar = 0, "T > Tf 
fits") could be well fit with the "stretched exponential" 
relaxation function 

P{t) [= P,{t) for T<Tf]= cxp[-(A<)^] , K <1, (3) 

often used^^ to model an inhomogeneous distribu- 
tion p{W) of exponential rates WM- The relaxation 
time 1/ A gives the time scale of the relaxation and the 
power K controls its shape: broader distributions of re- 
laxation rates are modeled by smaller values of K. For all 
fits K was found to be significantly smaller than 1, indi- 
cating that the spin fluctuations are inhomogeneous over 
the entire temperature range. At low temperatures and 
low fields the amplitude As of the slowly-relaxing compo- 
nent is approximately A/3, as expected from quasistatic 
LF-//SR relaxation in a powder sample.Ji 

Figure [2] gives the temperature dependencies of the pa- 
rameters from fits to data from both single-crystal and 
polycrystalline samples. For the two samples w^(T) and 
K reproduce well. With increasing temperature de- 
creases and vanishes at Tf = 8.5 ± 0.5 K, in agreement 
with the kink in the bulk susceptibility.^ A mean- field 
temperature dependence [Brillouin function for 5" = 1, 
curve in Fig. [H^a)] describes the data to within experi- 
mental uncertainty. 
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FIG. 2: (color online) Temperature dependencies of fit param- 
eters in NiGa2S4, ^loHl ~ 2 mT. Filled symbols: polycrys- 
talline sample. Open symbols; single-crystal sample, H_t | 
[001]. Red symbols: (a): oscillation frequency a)^/27r. (b): 
static relaxation rate Xr- Blue symbols: T <Tf fits (see text); 
(b) dynamic rate As, (c) stretching power K. Green symbols: 
T > Tf fits; (b) A, (c) K. Inset to (b): T > Tf, A/T^^^ vs 
l/T. Line: A/T^/^ oc exp(T7r) (Ref.[i3), T" = 114(7) K. 



The dynamic relaxation rates in the polycrystalline 
sample are higher below Tf and lower above Tf than in 
the single crystal but follow the same trends [Fig. ^b)]. 
As T approaches Tf from above A increases strongly, 
suggestive of critical slowing down of paramagnetic-state 
spin fluctuations. The data can be fit to A/r3/2 cx 
exp(r7T) [inset of Fig. [DJb)] with T* 114(7) K. This 
form is expected for critical fluctuations of Heisenberg 
spins on a 2D lattice^^ with T* — 27rp<,, where for 5 = 1 
the spin-wave stiffness ps « J. From 9w = zJS{S +l)/3 
with z = 6 Ni nearest neighbors, we estimate T* = 
27rps « 125 K, in good agreement with the measured 
value. The differences in A and K ior T < Tf and T > Tf 
fits near Tf [Figs. (H^b) and (c)] probably reflect a distri- 
bution of freezing temperatures. 

With decreasing temperature below Tf Ag decreases 
slightly, but for both samples remains > 1 /^s~^ down to 
2.1 K « Ty/4. This strong relaxation is consistent with 
the loss of ^^Ga NQR signal in the same temperature 
range 
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FIG. 3: (color online) Representative (a) early- and (b) late- 
time LF-/iSR asymmetry plots in single-crystal NiGa2S4 at 
T = 2.3 K and applied longitudinal fields Hl || [001]. 



Somewhat similar relaxation behavior has been 



observed in the triangular AFM NaCr02 (Ref. [IJ) with, 
however, no indication of critical slowing down above 
Tf 40 K. 

LF-/iSR asymmetry data from the single-crystal sam- 
ple at 2.3 K are shown in Fig. [3] for p.qHl in the 
range 1.13 mT-1.25 T. The rapid early-time relax- 
ation is decoupled (the late-time amplitude increases) 
for (J-oHl (Bloc). This decoupling is independent ev- 
idence for a distribution of quasistatic fieldsfii and is 
in contrast to the "undecouplable Gaussian" relaxation 
observed in /zSR experiments on SCGO (Ref. 19.). Fig- 
ured shows the field dependencies of A^ and K at 2.1 K 
(polycrystal) and 2.3 K (single crystal). The stretching 
power K x 0.4 is more or less independent of field. It 
can be seen that Ag is strongly suppressed by applied field 
fioHL ^ 10 mT -C (-Bloc). This suppression is not due to 
an effect of Hl on the niuon Zeeman splittingfi^ because 
the resultant niuon field /ioH^, -I- (Bioc) is still approx- 
imately (Bloc). Thus Hl must affect the Ni spin fluc- 
tuation spectrum directly. This contrasts sharply with 
the field-independent specific heat, which is unaffected 
to within a few percent by a field of 7 T. 

We conclude that NiGa2S4, which should be a good 
example of a 2D triangular Heisenberg = 1 magnet, 
is neither a "conventional" AFM, for which spin freezing 
would be accompanied by a decrease of the muon relax- 
ation at low temperatures, nor a singlet spin liquid, for 
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FIG. 4: (color online) Dependencies of (a) slow relaxation 
rate As and (b) stretching power K [Eq. (|3])] on applied longi- 
tudinal field Hl in NiGa2S4. Filled symbols: poly crystalline 
sample, T = 2.1 K. Open symbols; single-crystal sample, 
r = 2.3 K, Hl =|[ [OOf]. 



which there would be no spin freezing. The LF-/iSR data 
are consistent with a well-defined 2D phase transition at 
Tf « 8.5 K leading to a disordered AFM ground state, 
as found from previous ^^Ga NQR and fiSR studies4i^ 
This is in contrast to SCGO, which exhibits a field- 
independent specific heai*^ but no spin freezing The 
low-temperature relaxation rate As gives a lower bound 
on the frozen-spin correlation time Tc of 1/As ~ 10^^ s at 
2 K, more than three orders of magnitude longer than the 
lower bound from neutron scattering and seven orders of 
magnitude longer than the inverse exchange frequency. 
Thus the Ni spins are very viscous indeed, if not com- 
pletely frozen. The transition at Tf is unlikely to be due 
either to 3D coupling or to Ni-spin anisotropy, both of 
which are weak;— a topological transition associated with 
vortex binding^ is one candidate mechanism. NiGa2S4 



is chemically and structurally well ordered, raising the 
question of why the /xSR and ^^Ga NQR relaxation^ is 
spatially inhomogeneous. 

The persistent muon spin relaxation at low tempera- 
tures [Fig. [DJb)] is evidence for significant Ni spin fluc- 
tuation noise power at low frequencies. This is often 
observed in frustrated magnets without spin freezing^S 
but less frequently in spin-frozen stateSfi^i^iiS^ It in- 
dicates a high density of spin excitations that coex- 
ists with spin freezing at low temperatures A mean- 
field-like order parameter is also rare in 2D frustrated 
antiferromagnets^ and to our knowledge has never been 
seen together with persistent low-temperature dynamics. 

Easily the most remarkable result of this work is the 
contrast between the suppression of the low-temperature 
muon spin relaxation by Hl ^ 10 mT [Fig. |ll[a)] and the 
negligible magnetic field dependence of the specific heat 
up to fields more than two orders of magnitude greater.— 
There is no fundamental paradox here, since the entropy 
of a spin system is not explicitly dependent on the time 
scale of the fluctuations. But it is hard to understand 
how strong muon spin relaxation, usually associated with 
most or all of the spin excitation degrees of freedom of 
a system, could change so drastically with field without 
the slightest signature in the entropy. Reconciling the 
LF-/iSR and specific-heat results is clearly an important 
task for theory. 
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